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Abstract 

A new class of exact solutions to the axisymmetric and stationary vac- 
uum Einstein equations containing n arbitrary complex parameters and 
one arbitrary real solution of the axisymmetric three-dimensional Laplace 
equation is presented. The solutions are related to Jacobi's inversion prob- 
lem for hyperelliptic Abelian integrals. 
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1 Introduction 



The inversion problem for hyperelliptic integrals (integrals over a rational function 
R[x, W(x)] with W(x) being the square root of a polynomial of order 2n + 1 or 
2n + 2, n = 2,3,4,...) was formulated by Jacobi in 1832 [1]. Its solution in 
the ultraelliptic case (n = 2) was presented by Gopel [2] and Rosenhain [3]. 
They generalized Jacobi's theta functions to theta functions depending on two 
variables. Riemann and Weierstrass solved the general problem in terms of theta 
functions of n variables, cf. [4], [5], [6]. 

In the course of solving a physical problem within Einstein's theory of general 
relativity [7] , [8] — the problem of a rigidly rotating disk of dust — we were lead 
exactly to the ultraelliptic case of Jacobi's inversion problem [9]. In this letter 
we will show that a whole class of solutions to the axisymmetric and stationary 
vacuum Einstein equations may be associated with Jacobi's inversion problem in 
the general case. 



As is well known the axisymmetric and stationary vacuum Einstein equations 
may be reduced to the Ernst equation 



for the complex function f(p,() called the Ernst potential [10], [11]. The opera- 
tors A and V have their usual three-dimensional meaning applied to an axisym- 
metric function depending on the cylindrical coordinates p and £: 



2 The class of solutions 



(K/)A/ = (V/) 2 



(1) 



A 




(2) 
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Any solution of our class is given by 



f = exp < 



E 

m=l 



K n dK 



W 



- U r , 



(3) 



with 



W 



(K + iz)(K - iz)H(K - K 3 ){K - Kj). 
i=i 



(4) 



The K m (m — 1,2, ... ,n) are arbitrary complex constants; a bar denotes complex 
conjugation. The coordinates p and ( are combined to the complex coordinate 



z = p + iC 



(5) 



The upper integration limits (m— 1, 2, . . . , n) are functions of p and ( and 
have to be determined as the solution of the following inversion problem: 



?<"( m ) 
n - K^dK 



S / 

771=1 K 



W 



Uj, j = 0,1,2,..., n-l. 



(6) 



The quantities Uj (j — 0, 1, 2, . . . , n) in (6) and (3) are defined as follows: Let uo 
be any real solution to the (axisymmetric) Laplace equation 



Amq = 0. 



Then the Uj (j — 1,2, ... ,n) are real solutions to 



m i,z = ^tt- + zUj- liZ , j = l,2,...,n 



(7) 



(8) 



where an index ", z" denotes partial derivation with respect to z. As a conse- 
quence, all Uj are solutions to the Laplace equation. The calculation of the Uj, 
starting from u , is simply a (path-independent) line integration. If one starts 
from a solution u of the Laplace equation that vanishes at infinity one may 
choose the integration constants such that all Uj vanish at infinity as well. In 



this case one obtains / — > 1 at infinity, i. e. an asymptotically flat solution to 
the Einstein equations. It should be noted that the integrations in (3) and (6) 
have to be performed along the same curves in the two-sheeted Riemann surface 
(of genus n) associated with W(K) defined by (4). Using a result of Riemann 
[cf. [4], p. 311, Eq. (5)] it is possible to express / in terms of theta functions 
of n variables. On the symmetry axis, p — 0, the two zeros — iz and iz of the 
polynomial W 2 coincide. This leads, for n > 2, to expressions for the axis val- 
ues of the Ernst potential in terms of theta functions of n — 1 variables. Hence, 
the class of solutions considered cannot be obtained by applying usual Backlund 
transformations repeatedly to an arbitrary Weyl metric. 

The proof that / given by (3) - (8) satisfies indeed the Ernst equation (1) 
consists of two steps: 

(i) Combining Eqs. (8) and their complex conjugates with (6) and (3) one obtains 






(10) 



and 




(11) 




(12) 



where 




(13) 



Purely algebraic manipulations of Eqs. (9)-(12) lead to 



( ln /),* K (j) +iz II K (j) _ K{m ) ' 



(14) 



(to ^ j) 



4 



Kf ( ln f),z K(j) _■- II K (j)_ K {m) 

m = 1 

(m ^ j) 

The integrability conditions 



(15) 



yield, for all j, 



with 



K U)_ = K U) 

± ^-,zz ±l ,zz 



A(ln/) = a(Vln/) 2 



n 



m=l 



n 



(X(m) + iz p(m) _ ^ l = l i fl-(m) _ K (j) 



> . 



(16) 



(17) 



(18) 



(ii) Due to the reality of the Uj, it can easily be verified from Eqs. (6) that the 
and (m — 1, 2, . . . , n) are related such that they form the 2n zeros of 
a characteristic polynomial 



' n— 1 



U(K- K 3 )(K - Kj) - (K + \z){K -iz)\Y, <*i& 

j=l \j=0 

= (l-al_ 1 )f[(K-K^)(K-K^). 

3=1 



(19) 



The (purely imaginary) quantities a 3 - may be determined from the system of 
linear algebraic equations 



n—l 



, , w( m ) 
U (KW + iz)(Ki m )-iz) 



m — 1, 2, . . . , n 



(20) 



following from (19) for K = (m = l,2,...,n) and (4). In particular, it 

turns out that 

a n -i = a (21) 
with a defined according to (18). Using (19) together with (3) one finds 

f-f 



f + f 



OLn-l 



(22) 



and consequently, because of (17), 

A(ln/) = ^|(Vln/) 2 (23) 
which is equivalent to the Ernst equation (1). 

3 Discussion 

For n — the class of solutions (3) reduces to Weyl's class with / = exp(— it ). 
The case n — 1 leads to solutions in terms of elliptic functions. A nontrivial and 
physically meaningful example for n = 2 is the solution of the problem of a rigidly 
rotating disk of dust [9]. This example provides some justification for the hope 
that the class of solutions presented will lead to further physical applications in 
the context of the rotating body problem. 

An interesting mathematical question concerns the relation to Korotkin's 
finite-gap solutions [12]. It seems to us that those solutions can be obtained 
for a particular choice of the potential function w and form therefore a subclass 
of the solutions presented here. 
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